Spectral Properties of Directed Random Networks with Modular Structure 
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We study spectra of directed networks with inhibitory and excitatory couplings. We investigate in 
particular eigenvector localization properties of various model networks for different value of corre- 
lation among their entries. Spectra of random networks, with completely uncorrelated entries show 
a circular distribution with delocalized eigenvectors, where as networks with correlated entries have 
localized eigenvectors. In order to understand the origin of localization we track the spectra as a 
function of connection probability and directionality. As connections are made directed, eigenstates 
start occurring in complex conjugate pairs and the eigenvalue distribution combined with the local- 
ization measure shows a rich pattern. Moreover, for a very well distinguished community structure, 
the whole spectrum is localized except few eigenstates at boundary of the circular distribution. As 
the network deviates from the community structure there is a sudden change in the localization 
property for a very small value of deformation from the perfect community structure. We search for 
this effect for the whole range of correlation strengths and for different community configurations. 
Furthermore, we investigate spectral properties of a metabolic network of zebrafish, and compare 
them with those of the model networks. 
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I. INTRODUCTION 

Statistical properties of large random matrices have 
been investigated intensively since the early 1950s and 
have turned out to be valuable tools for both qualita- 
tive and quantitative descriptions of spectral properties 
of complex systems Examples include, complicated 
interactions in nuclei, quantum mechanical aspects of 
chaos, and the stock market. Recent success has been 
achieved in systems having an underlying network struc- 
ture [31 . In applications of random matrix theory to phys- 
ical problems, it is generally assumed that the details of 
the physical systems are less important for many statis- 
tical properties of interest. Often it turns out that im- 
portant statistical properties such as the distribution of 
eigenvalues or the spacings of energy levels in quantum 
systems are well described by the respective properties 
of random matrices that respect the same symmetries as 
the physical system [1] . 

Recent investigations of statistical properties of eigen- 
values of networks show that bulk eigenvalue spacings 
follow Gaussian orthogonal ensemble (GOE) predictions 
of random matrix theory (RMT) d-Q, whereas the eigen- 
value distribution itself depend on network properties 
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[5|-|9|. These analyses show that the spectral distri- 
butions of bulk eigenvalues of random networks follow 
Wigner's semi circular distribution, for scalefree networks 
the bulk eigenvalues follow a triangular distribu- 
tion, for smallworld networks (Tl| they form a multi peak 
distribution, and for realworld networks they show one 
of the above mentioned distributions with a high peak 
at zero eigenvalues Q . Further investigations of the 
community structure in networks show the transition to 
Wigner's semi circular distribution as the community 
structure is destroyed. Even a very small deformation 
from the community structure gets captured by the eigen- 
values distribution [l2}]. 

All the above investigations have been performed for 
undirected networks leading to a symmetric adjacency 
matrix and for which eigenvalues are real. Many com- 
plex systems are described by directed networks [1^, for 
example, the World Wide Web, neural networks, pro- 
tein interaction networks and many social networks are 
directed or asymmetrically weighted, and hence lead to 
complex eigenvalues of the corresponding adjacency ma- 
trices. One of the most important features shown by real- 
world networks is the existence of modular or community 
structures The study of community structures 

helps to elucidate the organization of networks, and even- 
tually could be related to the functionality of groups of 
nodes. Regardless of the type of realworld networks in 
terms of the degree and other structural properties [11] , 
it is possible to distinguish communities throughout the 
whole networks [l^. However, the question of the defini- 
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tion of the community is problematic. Usually commu- 
nity is assigned to the nodes that are connected densely 
among themselves, and are only sparsely connected with 
other nodes outside the community. 

So far most of the theoretical investigations and ap- 
plications of RMT have focused on symmetric matrices. 
Asymmetric matrices or non-Hermitian matrices are less 
well understood. A classic result for non-Hermitian ma- 
trices in random matrix theory is Girko's circle law [l7j 
which states that, for large N, the eigenvalues of an 
N X N asymmetric matrix lie uniformly within the unit 
circle in the complex plane, if the elements are chosen 
from a distribution with zero mean and variance 1/N. 
When partial symmetry is included, the circle changes 
to an ellipse [l8|. Recently, numerical studies on spec- 
tral properties of asymmetric matrices have been done 
in the context of Google networks [Ull^. Additionally 
the spectral density of non-Hermitian sparse matrices has 
been analyzed in the Ref. ^21i |. 

In this paper we evaluate the applicability of RMT for 
the distribution of eigenvalues of directed networks where 
entries in the corresponding matrix take values motivated 
by inhibitory and excitatory coupling between nodes. In 
particular, we investigate spectral properties of directed 
networks having a community structure. The matrix cor- 
responding to these networks is different from the non- 
Hermitian random matrices studied under a RMT frame- 
work. The main differences are that (a) the matrix is not 
random^ i.e. the position of entries in the matrix depends 
on the underlying network structure; (b) real world net- 
works and model networks are sparse and consequently 
most of the entries in the corresponding matrix are zero; 
and (c) most importantly, the entries in the matrix take 
values 0, —1 and 1. These entries are motivated by the 
excitatory and inhibitory synapse between neurons lead- 
ing to positive and negative values of the connections 
between them. 

The spectra of matrices having entries 0, 1 and — 1, 
known as Seidel matrices, have been investigated inten- 
sively in the literature. In graph theory, the Seidel adja- 
cency matrix of a simple graph G is a symmetric matrix 
having at the diagonal entries, Gij = — 1 if node i is 
neighbor of node j, and other entries of this matrix tak- 
ing the value 1. For Seidel matrices Gij = Gji. The 
multiset of eigenvalues of this matrix is called the Seidel 
spectrum. The Seidel matrix was introduced by van Lint 
and Seidel [12] and extensively exploited by Seidel and 
co-workers. Seidel matrix provides insight to the proper- 
ties of graphs, thus leading to several applications such 
as modeling social situations, spin glasses in physics, and 
clustering of data [2^ . 

Note that, except for the matrix entries being binary 
1,-1 and 0, Seidel matrices are completely different from 
the matrix representation considered in the present pa- 
per. The first difference between Seidel matrices and the 
matrices investigated herein is that the latter are directed 
(asymmetric) matrices, i.e., Aij ^ Aji depending upon 
the value of r, a measure of directionality or correlation. 



Consequently, we get complex eigenvalues following the 
universal Girko's law for the case when the matrix is com- 
pletely uncorrelated. The second most important differ- 
ence is that we have Aij = for the nodes that are not 
connected, i.e., vertices that are not neighbors, whereas 
Seidel matrices have the entries 1 for non neighboring 
nodes. This second property of our matrices again leads 
to a very crucial difference called 'sparseness', which af- 
fects the spectral properties to a great extent in compari- 
son to the matrices corresponding to fully connected net- 
works. The third important difference between the ma- 
trix we have considered in this paper and the Seidel ma- 
trices we found in the literature is the underlying struc- 
ture: Seidel matrices are mostly investigated for regular 
graphs (such as strong regular graphs and two-graphs) 
[2J|, whereas the matrices that we are investigating in 
the paper are random with some structure, correspond- 
ing to ER graphs, BA scalefree graphs, modular graphs 
and modular realworld networks. There are cases in the 
literature where the Seidel matrix of ER graphs and re- 
alworld graphs have been investigated; in the following 
we provide details of these studies only, as they seem to 
be close to our studies at least in terms of the structure 
of the underlying network ( Seidel matrices of random 
graphs have been studied mainly for the graph bisection 
problem.) Reference [25!| shows through numerical anal- 
ysis that Seidel representations produce the smallest cuts 
for random ER graphs than other representations such as 
the adjacency and Laplacian matrices. 

Spectral densities of directed random networks with 
entries 0, 1 and —1 have also been studied analytically 
and numerically in Ref. [1^. In Ref. [1^ eigenvalues 
spectra of large random matrices with excitatory and 
inhibitory columns drawn from a distribution with dif- 
ferent means and equal or different variances have been 
computed. Interestingly, the spectra of these random 
matrices (with and 1, —1 entries) also follow Girko's 
circle law. 

The main goal of the present work is to investigate the 
changes in the spectral properties of the corresponding 
matrix as the underlying directed network deviates from 
the perfect community structure. Furthermore, in or- 
der to understand the origin of delocalized eigenvectors, 
we track the spectra as directionality is introduced into 
the networks, i.e. connections are made directed with 
probability t. Apart from the ispectral distribution we 
study eigenvectors localization properties of model net- 
works and a realworld networks. 



II. MODEL 

We generate networks with a community structure as 
follows. First, m Erdos-Renyi random networks [l^] with 
connection probability p are constructed; the spectral be- 
havior of the matrix corresponding to each of these sub- 
networks (blocks) separately follows GOE statistics [l^ . 
The matrix corresponding to the full network would be 
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an m block diagonal matrix. Then, with probability pr 
connections are rewired as follows. A pair of nodes from 
a sub-network is chosen randomly. Subnetwork is also 
chosen with equal probability. The connection between 
these two nodes is deleted, and a new connection between 
one of the node in the chosen pair and a randomly chosen 
node outside to this sub-network is introduced. This con- 
figuration leads to an m block matrix with, blocks having 
entries 1 with probability p, and off-diagonal blocks hav- 
ing the entries 1 with probability p^- The above networks 
can be casted in the following form: 

A^Aq + Ap^ (1) 

where Aq is the direct sum of all the blocks {Aq = 
^1 ® ^2 © ■ • • Am)- The directionality of the edges is in- 
troduced as follows. A node is assigned equal probability 
of being inhibitory or excitatory. All edges starting from 
an excitatory node would take value 1, and those starting 
from an inhibitory node would take the value —1. This 
means that the corresponding adjacency matrix A has 
entries as follows: if the ith node is excitatory, Aij = 1 
for all the j nodes that have connections with i, and if 
the ith node is inhibitory, Aij = — 1. The corresponding 
matrix has following two properties: (a) \Aij\ — \Aji\, 
and (b) a row has all entries either 1 or —1. This con- 
figuration would lead to zero mean ^ Aij ^ 0), and 
standard deviation ^ r. We denote the eigenvalues 
of network by A„, n = 1, . . . , miV, where N is the size of 
sub-network, and m is the number of sub-networks. Note 
that the size of each sub- network may be different, but 
for simplicity we consider equal size here. 

Connections are rewired with the probability < < 
1, where Pr defines the ratio of inter-community (Ninter) 
to intra-community {Nintra) connections. Starting with 
two ER sub-networks of dimension N, a fraction pr of 
the connections is rewired. The average degree of net- 
work for each pr remains the same as that of the initial 
configuration. Nodes are inhibitory or excitatory with 
equal probability, which implies an approximately equal 
number of 1 and —1 entries randomly distributed in the 
corresponding matrix. For pr ^ 0.5, the whole struc- 
ture becomes a directed random network, and for ~ 1 
which indicates that all the connections are rewired, the 
network takes the form where nodes in ^1(^42) are con- 
nected to those in ^2(^1). 



we use the inverse participation ratio (IPR) defined by. 



where '>pi{j),j — 1,...,A^ are the components of the 
eigenvector ip^. The meaning of / is illustrated by two 
limiting cases : (i) li = N for a vector with identical 
components ipi{j) = whereas (ii) li — 1 for a vec- 

tor with one component V'i(j) = 1 ^-nd the remainders 
zero. Thus, the IPR gives the effective number of nodes 
on which an eigenstate is localized. 

Figure ^ plots eigenvalues together with the IPR 
for the model described in the section [TTl Subfigures 
correspond to different rewiring probabilities. For each 
rewiring probability eigenvalues are distributed homoge- 
neously in a circular region. Gray shading (color) de- 
notes the value of IPR for the corresponding eigenstate. 
For Figure (P a) with ~ 0, which corresponds to the 
network having two communities connected with one sin- 
gle connection, large number of eigenstates have approxi- 
mately equal IPR value close to the bottom of the shading 
(coloring). Few eigenstate lying at the boundary of the 
circle and a few near the real axis have low IPR values, 
indicating high participation in the localization. 

Figure ([T]) plots the results upto rewiring probability 
p = 0.01, as rewiring is increased further the spectral dis- 
tribution remains exactly the same as for p = 0.01 case 
except that the light gray (green) hazy circular part be- 
comes bigger and denser indicating that more eigenvalues 
are delocalized. As rewiring keeps on increasing, leading 
to the destruction of the community structure more and 
more eigenvalues join the light gray (green) big circle, and 
for Pr = 0.1 except for a few eigenvalues lying at the end 
of the real and imaginary axes all are delocalized. Ad- 
ditionally, the eigenstate near the real axis has relatively 
lower IPR values. One eigenvalue lying at the extreme 
end of the real axis is the most localized. For a further 
increase in rewiring probability there is not much of a 
change in the localization properties of eigenstates. The 
value Pr — 0.5 corresponds to the same expected number 
of intra and inter-community connections. At this value 
the network is a complete random network. A further 
increase in pr takes the network toward a bipartite type 
of structure. 



III. SPECTRA 

In the following we present the results for N — 1000, 
m = 2, and p = 0.02. The average degree of the network 
can be calculated as < fc >~ pN = 20. The eigenvalues 
spectra of Ai and A2 forp^ = follow Girko's law of 
non-Hermitian matrices |17| . We numerically diagonal- 
ize the adjacency matrix of the network to obtain a set of 
eigenvalues = Rk + Hk, and eigenvector ipi. To char- 
acterize the localization properties of the eigenvectors ipi, 



A. Fine rewiring 

Figure ([T|) shows that the localization properties of a 
network change completely even for a very small value 
of pr, as small as Pr = 0.1. After this value no visi- 
ble changes are seen in the localization properties of the 
bulk of the spectra. To gain insight into the localiza- 
tion to delocalization transition we track the spectra for 
each rewiring between the communities (Eq. [T|). First 
four subfigures (Figure [T] a - [T] d) plot the spectra along 
with the localization value for a very small number of 



4 



5 Cay. 



5 Ch). 



-Si CcX 



— 5 O S 



O S —5 O 5 



rcw=10 



f 

Cd) 



—5 O 5 




5 O 5 



FIG. 1: (Color online) Spectra with IPR for networks start- 
ing with two communities of size TV — 1000 and then rewiring 
connections between communities with different probabilities. 
Average degree is < fc 20. The real and imaginary parts of 
the eigenvalues are plotted on the x and y axes, respectively, 
(b) One and (c) five rewirings are plotted for an ensemble 
average of 100. For other rewiring probabilities the spectra 
are plotted for an ensemble average of 10. The rewiring prob- 
ability is denoted by pr- The value of the IPR varies from 
(a)-(f) / = [ dark gray, blue, patches at the corners] to (b)- 
(f) / = 1100 [dark gray, red, points inside the circular region]. 
Light gray regions inside the circle in (a) - (c) depict less de- 
localization than the dark gray patches at the corners. Light 
gray area inside the circular regions and the horizontal line at 
X axis of (d)- (e) are even less delocalized than the peripheral 
light gray areas that have the same level of delocalization as 
the interior of (a)-(c). 



rewiring All figures are plotted for an ensemble over 100 
random initial conditions for the networks. Figure ([T] a) 
plots the spectra and IPR value for a perfect community 
network without any connections between them. This 
configuration leads to two 

disconnected components. With only one rewiring be- 
tween two communities, we already start observing the 
delocalization of few eigenstates (Figure [T] b). Each ad- 
ditional rewiring between communities leads to the delo- 
calization of more eigenstates. For as few as 10 rewiring 
to the initial configuration (Figure [T] a) a large number 
of eigenstates are delocalized (Figure [T] d) . One can no- 
tice in Figure ((J) that delocalized eigenstates form a hazy 
ring around center [2^. Eigenstates lying at the periph- 
eral of the circle and those corresponding to the real axis 
delocalized in last for very large value of rewiring prob- 
ability between the communities. For pr = 0.01, eigen- 
states near real axis and with large absolute values are 
those which still maintain localization. As value of Pr 
increases further except few eigenstates having largest 
absolute values lying at real and imaginary axis are de- 
localized. Bulk middle part delocalized completely. As 
value of Pr increases further, bulk middle part remains 
same, only eigenstates with largest absolute eigenvalues 
at imaginary axis are more delocalized. For pr ^ 0.9, the 
spectra remains similar to that of value pr — 0.1. 



We would like to note that spectral properties of m 
coupled GOE random symmetric matrices have been in- 
vestigated in details. They have been used to understand 
the behavior of quantum graphs (28j . 



IV. COMPARISON WITH RANDOM 
NETWORKS 

In order to gain insight into community induced spec- 
tral changes, we present results of random networks. In 
the following we first provide the eigenstate properties of 
random networks having different average degrees, and 
then track the spectral properties as parameters of r. 



A. Completely uncorrelated (directed) network 

First we consider completely uncorrelated random net- 
works, for which r = and the connection probability is 
p. Since there is equal probability of a node being in- 
hibitory or excitatory, the expected number of 1 and — 1 
entries is approximately same. The mean and variance 
of this network can be calculated as /i = and cr^ = p, 
respectively. 

Figure ^ plots eigenvalues together with the IPR for 
random networks with different connections probabili- 
ties. 

As we do not consider self connection, the diagonal 
elements of the matrix are zero, and hence eigenvalues 
are distributed around origin 0. Eigenvalues are dis- 
tributed in circular region of radius Np{l — p). For 
larger values of p, few pairs of eigenvalues get separated 
from the circular bulk. The bulk circular region still lie 
between — Np{\ — p), ^ Np{l — p). For p = 1 which 
corresponds to globally coupled network, there is one 
conjugate pair of eigenvalues with absolute value -s/TV 
which are nondegenerate, rest of eigenvalues are degen- 
erate with values 1 and —1. 

Different connections probabilities lead to different ex- 
pected average degree of the network. For low connection 
probability, p = 0.005, one gets localized eigenstates at 
four corners of the real and imaginary axis. Apart this 
four corners the eigenstates on real axis are also more lo- 
calized. As connection probability increases overall dis- 
tribution of eigenvalues remains same, i.e. homogeneous, 
except for larger p few pairs of eigenvalues get separated 
from the bulk circular region [l^ . Except eigenstates ly- 
ing on real axis, other eigenstates are toward dark gray 
(blue) denoting delocalization. The isolated eigenstates 
are maximal delocalized, which implies that eigenstates 
having large absolute eigenvalues are more random than 
the bulk part. Very small values of p yields a delocalized 
spectra (Figure [2] a), which may be due to the sparseness 
of connections. 

Here we would like to mention that, the above behav- 
ior is for the special arrangement of 1 and —1 entries 
in the matrix, a row has all 1 or —1 entries depending 
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FIG. 3: (Color online) The IPR for (a) random network with 
the community as in Figure [T] as a function of the connection 
probability between communities pr, N = 2000 and average 
degree < k >= 20 (b) a random network, A*' = 2000 as in 
Figure ([2]), as a function of p or average degree < k >. 



FIG. 2: (Color online) Spectra with the IPR for directed ran- 
dom networks having different connection probabilities p. The 
size of the networks A'^ — 2000, and spectra are plotted for 
10 random realizations for each value of p. The real and 
imaginary parts of the eigenvalues are plotted on the x and y 
axes, and gray shading (color) denotes the value of the corre- 
sponding IPR. The value of the IPR varies from / = [dark 
gray, blue, patches in corners of (a)] to / = 2000 [scattered 
dark gray, red, dots in (d)]. Light gray areas inside the circu- 
lar regions in (a)-(d) are less delocalized than the dark gray 
patches in (a). Light gray dots outside the circular region in 

(c) are not as delocalized as the light gray regions described 
earlier. Darker gray dots outside the circular region in (c) and 

(d) are not as delocalized as the light gray dots in the same 
region and are more delocalized than the dark gray dots in 
(d), which are most the localized. Note the scale of the gray 
shading (color) axis for this case, which is different from other 
figures, (a)-(d) The horizontal line at x-axis shows the same 
level of localization as the value of delocalization close to the 
most delocalized eigenstates. 



upon whether the corresponding node is excitatory or in- 
hibitory. For a random network with entries 1 and — 1 
randomly distributed, the radius of circular region scales 
with average degree of the network, i.e. \^pN, and all 
eigenvalues lie within the bulk regions even for larger 
value of p including p = 1. 

For random symmetric networks (r = 1, in EqU]) with 
mean p and variance p{l — p), the distribution of the 
largest eigenvalue is given by the results in j29j . which 
shows that for constant p the distribution of the princi- 
pal eigenvalue of the random network is asymptotically 
normal with mean Np+ (1—p) and variance 2p{l — p) as 

— > oo . For asymmetric random networks the corre- 
sponding result for the distribution of largest eigenvalue 
has yet to be found, but in (30j it was shown that asym- 
metric network of and 1 entries and having mean p and 
variance p{\—p), the principle eigenvalue asymptotically 
converges to Np. For matrices with entries 1 and —1, the 
mean and variance differ from above, as mean and vari- 
ance for completely uncorrelated random network (r = 0) 
with parameters N, p would be and p respectively. 

Though the IPR is, in general, used to understand the 
localization property of individual eigenstates, the sum 



of all li 

<I>^<Y,h>, (3) 

i 

where bracket means ensemble average, can be used as a 
measure for localization properties of whole spectra alto- 
gether [28| . We calculate the total value of IPR < / > in 
order to get a measure of localization for the whole net- 
work. Figure ([3]) plots < I > for a random network with 
a community structure for an increasing value of the con- 
nection probability between communities. In a straight 
forward manner, the figure shows delocalization in the 
whole spectrum with deformation from the community 
structure. The value of < / > increases, which means 
that overall the spectrum becomes more delocalized un- 
till a certain value of rewiring probability pr ^ 0.005 and 
then becomes saturated, as shown in Figure (jlj. The 
localization properties the of spectra are barely distin- 
guishable from the community structure for larger value 
of deformation. The figure is plotted for network size 
N — 2000 and average degree < k >= 20. We also 
plot the IPR ([3]) of random networks for various value 
of p or an average degree < k >. the network size al- 
ways remains the same, N = 2000. The range for p is 
from p — 0.001 to 0.9. We do not plot the graph below 
p ^ 0.001, since for N — 2000, smaller values of p yield 
several disconnected clusters. N = 2000 and p ~ 0.001 
yield an average degree < k >^ 20. 



B. Tracking from correlated (symmetric) to a 
completely uncorrelated (asymmetric or directed) 
network 

In this section we investigate the origin of localization 
or delocalization of eigenstates by tracking the spectra as 
connections are made undirected. We start with a net- 
work having all nodes excitatory leading to a symmetric 
network, now with probability (1 - t)/2 some nodes are 
made inhibitory, r = corresponds to the case when 
half of nodes are made inhibitory, leading to a completely 
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FIG. 4: (Color online) Spectra with the IPR for random net- 
works having different values of r. (a)-(f) The size and av- 
erage degree of the network are A'^ = 2000 and < k >= 10 
respectively. For a symmetric network, t — 1, and r = 
corresponds to complete asymmetric network (see the text). 
The value of the IPR varies from / = [dark gray, blue, at 
the bottom of shading scale] to / = 1300 [dark gray, red, at 
the top of the scale] . Dark gray points at both ends of the line 
at a; = and dots at the periphery of the two circular region 
on both sides of the real axis correspond to 7 = 0, whereas 
a dark gray dot separated by other eigenvalues in (a) and 
(b), and dark gray patches inside the circular region in (c)-(f) 
correspond to 7 = 1300. Gray shading in (b) denotes differ- 
ent levels of delocalization from 7 = (dark gray) to 7 = 600 
(light gray) , whereas in (c) and (d) , gray shading in the circu- 
lar region denotes the localization value from I = 1000 (dark 
gray) to 7 = 600 (light gray), with other gray shading for the 
intermediate values. Light gray regions in the middle of the 
real axis and toward the top and bottom of the circular re- 
gions in (c)-(f) correspond to the middle of the scale, whereas 
the darker gray region in (c) and (d) correspond to the value 
7 = 800 in the gray scale. 

uncorrelated network. As r fraction of connections are 
made symmetric, the mean and variance take form as: 

T = 1 corresponds to a symmetric network which we call 
as a complete deviation from the directed network, i.e. 
Aij = Aj^i Mi&ij. We call "deviation from direction- 
ality" because as one deviates from completely directed 
networks, spectra deviate from the circular structure. For 
a network, where nodes are inhibitory or excitatory with 
equal probability, < Aij >— and < A^j >^ p. 

Figure (j4]) plots the spectra along with the IPR values 
for random networks for different values of r. The largest 
eigenvalue lies on real axis well separated from the bulk. 
Isolation of largest eigenvalue from rest is well known ob- 
servation for real symmetric random matrices having non 
zero mean fsH .It was proved in [sij that such an isolated 
eigenvalue exists for a large matrix if the mean value of 



the elements is not zero and rather carries substantial 
fraction of the root mean square of the elements. Fur- 
ther more, the associated eigenvector would have all its 
components close to equality. 

Now let us first consider the case with r = 1. The 
corresponding matrix would be a symmetric matrix with 
pN{N— 1) entries being 1 and rest being 0, the mean and 
variance of this matrix can be calculated as fj, = p and 
cr^ — p{l—p) respectively (Eq.|4]). The largest eigenvalue 
of the matrix scales as pN, which is equal to the aver- 
age degree of the network. Rest of the eigenvalues are 
homogeneously distributed in a circular region of radius 
pN{l —p). As connected are made directed, the mean ji 
decreases, and for r = 0, mean takes value 0, and vari- 
ance is simply given by p which is connection probabil- 
ity. As connections are made directed, eigenvalues start 
appearing in complex conjugate pairs, and we can di- 
vide spectra into two parts: (A) part of spectra with real 
eigenvalues, and (B) the part of spectra having complex 
conjugate pairs of eigenvalues. For r — 0.8, which means 
20% of connections are made directed, many eigenvalues 
are still in group (A), i.e., they lie on real axis. Complex 
conjugate pairs of eigenvalues form an oval shape that is 
well separated from real axis (see Figure (|4]). The local- 
ization properties of eigenstates from both the groups do 
not change much, as bulk part which is separated from 
the real axis is still delocalized, only middle part of eigen- 
states lying at the real axis (group (A)) delocalized than 
a completely symmetric network. As more connections 
are made directed, complex conjugate pairs of eigenval- 
ues start appearing with delocalized states. Bulk of the 
complex conjugates pairs lying at both sides of the real 
axis and forming oval shape (group (B)) are delocalized, 
only few eigenvalues in this group with largest absolute 
value remain localized. Eigenvalues forming group (A) 
also show similar features as group (B), eigenstates with 
larger eigenvalues at both ends remain localized, whereas 
rest of the eigenstates are more delocalized than earlier. 

As more connections are made directed (value of r is 
decreased), size of group (B) keeps increasing, and conse- 
quently two parts of eigenstates come closer. For r = 0.2, 
group (B) is no longer separated from the eigenstates ly- 
ing on the real axis (group (A)). 

Figure ^ plots the total < I > (Eq[3]) for a partic- 
ular value of connection probability p = 0.01 and net- 
work size N — 2000. Starting from a symmetric network 
(t = 1), connection are made directed with probabil- 
ity (1 — r). Overall spectral delocalization increases as 
connections are made directed leading to occurrence of 
complex conjugate pairs of eigenvalues in spectra. De- 
localization reaches its maximum for t ~ value which 
corresponds to a complete uncorrelated or directed net- 
work. For networks having only positive entries, there is 
not much changes in the localization properties till value 
T < 0.8, after this there is a sharp decay in < J > indicat- 
ing sudden localization transition. For values between 
and 0.8, spectral distribution smoothly changes from cir- 
cular to elliptical, localization property does not show 
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FIG. 5: (Color online) The total IPR as a function of r. The 
data are plotted for TV = 2000, < >= 20 and for 10 en- 
semble averages of the random networks, (a) Networks with 
inhibitory and excitatory nodes leading to 1 and —1 entries in 
the corresponding matrix, and (b) random network with only 
excitatory nodes leading to only 1 entries in the corresponding 
matrix. 



any rich structure, as shown for the previous case. 



C. Random networks with positive entries 

In the above we investigated networks with both 
excitatory and inhibitory connections. This section 
considers all nodes being excitatory, and consequently 
the adjacency matrix has all positive entries only. For 
asymmetric random network of and 1 entries, with 
mean p and variance p{l — p), the principle eigenvalue 
asymptotically converges to Np [30|- The rest of the 
eigenvalues are distributed in circular region with radius 
given by Np{l — p). In this section we present the 
spectral changes as network goes from an asymmetric 
(directed) corresponding to r = 0, to a symmetric 
network corresponding to r = 1. The mean and variance 
for networks with entries 1 would not depend on the 
value of T, and depends only on connection probability 
p as 



fj,=p 



(5) 



Figure ^ shows spectra along with localization prop- 
erty for different values of r. Spectral distribution 
smoothly goes from circular distribution to elliptical as 
value of T increases. For all values of r, the largest eigen- 
value is completely delocalized, and situated very well 
separated from the rest bulk part with value Xmax ^ Np, 
as for the case r = 0. Note that for each value of r, the 
expected number of entries 1 remains same for a fixed p 
value. 



V. NETWORKS WITH SCALEFREE 
COMMUNITY STRUCTURE 

In this section we consider that each community is 
modeled by a scalefree network instead of a random net- 
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FIG. 6: (Color online) Random network for various values 
of T but with only excitatory coupling. The network size 
— 2000 and the average degree < k >= 20. The value of 
the IPR varies from / = [light gray, blue, at the two ends of 
the spectra on the x axis] to / = 1200 [dark gray, red, circular 
regions] . the light gray region in the middle of x-axis in (f ) is 
not as delocalized as the light gray points toward both ends. 
The light gray horizontal line at x axis in (a)-(c), and the 
light gray patches at the corner of (b)-(e) have the same level 
of localization, which shows less delocalization than the light 
gray region in (f) described above. 



work. Each sub-network has size N — 1000 and average 
degree < k >— 20, and is constructed using BA pref- 
erential attachment rule [l^ . With probability Pr intra- 
community connections are rewired into inter-community 
connections. Connections are made directed in same 
manner as given in the section |TT1 

Spectra of the above construction are plotted in Fig- 
ure ([T]). The first difference between the spectra of these 
networks from those having a random community struc- 
ture (jni is that the eigenvalues are not homogeneously 
distributed in the circular region. The spectra are more 
dense around origin, and sparser near boundary. Also the 
value of largest absolute eigenvalue increases, as the high- 
est degree of the node in the current structure are very 
much larger than that of the previous structure (|lT]) . Lo- 
calization property overall remains same as those of the 
networks discussed in |lTl 



VI. COMMUNITY STRUCTURE : TRACKING 
FROM SYMMETRIC TO A COMPLETE 
ASYMMETRIC (DIRECTED) NETWORK 

In this section we investigate the origin of localization 
or delocalization properties of eigenstates for networks 
having distinguished community structure. Particularly 
we concentrate on the community dependent localization 
features as connections are made directed, and compare 
them with random networks features presented in the 
section IIVBI Starting with 171 — 2 random symmetric 
sub-networks with pr — 0.2 rewiring probability between 
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FIG. 7: (Color online) Eigenvalues with the IPR for networks 
having two scalefree sub-networks, and for different rewiring 
probabilities. The size of each community A'' = 1000 and 
the number of communities is m = 2. The average degree is 
< k >= 20. The rewiring probability is denoted by pr- The 
value of IPR varies from / = [dark gray, blue, scattered 
points at the peripheral of all sub-figures and a few points 
separate from the bulk in (d)-(f)] to 7 = 1200 [dark gray, red, 
patches in the middle of (b)-(f)]. Gray shading in (a) shows 
different levels of delocalization with dark gray being the most 
delocalized to lighter gray being less delocalized. Light gray 
points at the periphery in (b)-(f) are not as delocalized as 
the light gray region inside the circle of (a), whereas they are 
more delocalized than the light gray region inside the circle 
in (b)-(f). 



these sub-networks, connections are made directed with 
the probability r as explained in the previous section. 
Again, r = corresponds to a fully directed network 
such as explained in and t — 1 corresponds to the 
symmetric network. A symmetric random network with 
connection probability p has mean p and variance p{l—p) , 
whereas mean and variance of a fully directed network 
(made of 1 and —1) are and p respectively. Probability 
of directed connections (r) affects the mean and variance 
as given in Eq. |4l 

Figure ^ plots the eigenvalue spectra combined with 
the IPR as r varies from t = 1 (a perfect symmetric case) 
to T = (a uncorrelated case). Subfigure (a) corresponds 
to the configuration when almost 80% nodes are excita- 
tory leading to very large correlation among entries of 
the matrix. For very small number of directed connec- 
tions the spectra have mostly real eigenvalues, and eigen- 
states are localized. As number of directed connections 
increases spectra have more complex eigenvalues. Sub- 
figure (i) corresponds to the case where ~ 80% of nodes 
are inhibitory, which is exactly the same case in terms 
of correlations as shown for subfigure (a) except that the 
corresponding network has dominating —1 entries. 

As similar to the section IIVBI spectra can be divided 
into two parts, (A) part of spectra lying on the real axis, 
and (B) the part which form complex conjugate pairs. 
The appearance of complex eigenvalue seems to be an 
indicator of delocalization transition. For small number 
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FIG. 8: (Color online) Spectra with the IPR for networks 
having two random sub-networks, and for pr = 0.2 rewiring 
probability. The size of each community A'^ = 1000 and the 
average degree < k >= 20. The different probabilities of di- 
rected connections are plotted (see the text). Configurations 
are shown in which (a)-(d) excitatory nodes are dominant, 
(f)-(i) inhibitory nodes are dominant, and (e) there is equal 
probability of a node being excitatory or inhibitory. The value 
of the IPR varies from I — [dark gray, blue, patches at the 
periphery of both circular regions on two sides of the real 
axis in (a) and (i), and dark gray points toward both ends 
of the real x axis in (c)-(g)] to / = 1200 [dark gray, red, 
points well separated from bulk in (a)-(i)]. Gray shading in 
(a),(b), (h) and (i) denotes different levels of delocalization 
with dark gray (except two dots well separated from the bulk 
) denoting maximum and light gray showing minimum delo- 
calization. Light gray area inside the circular regions in (c) 
and (g) show less delocalization than the light gray area in 
(a) and (b). Dark gray area inside the circular regions in (c) 
- (g) show less delocalization than the similar dark gray area 
inside the circular region of (a). 



of directed connections (Figure [8] a) most of the eigen- 
values are real and hence localized. As more connections 
are made directed, complex pairs of eigenvalues start ap- 
pearing and eigenstates associated with these complex 
values start getting delocalized. For r = 0, which leads 
to a complete asymmetric or directed network, one gets 
complex pairs of eigenvalues covering the whole circle, 
and hence almost all eigenstates, except those associated 
with real eigenvalues and those with maximum absolute 
value, get delocalized. 

Now, as we compare the spectra of above construction 
with those of random networks, we observe one remark- 
able difference revealing two communities hidden in the 
network. For small number of directed connections (i.e. 
large values of correlation parameter t), there are ex- 
actly two isolated eigenvalues on the real axis well sepa- 
rated from the bulk. Figure (jS] a) shows that these iso- 
lated eigenstates corresponds to two largest eigenvalues 
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FIG. 9: (Color online) (a) Spectra of the zebrafish metabolic 
network, (b) Close-up of the spectra, (c) Spectra of the ran- 
dom network with same parameters (i.e. number of nodes, 
and number of directed and undirected edges) as zebrafish 
metabolic network. Gray shading (color) denotes the local- 
ization properties of the corresponding eigenvalues. Dark gray 
(blue) patches in (a)-(b) and dark gray dots at the end of the 
real axis as well at the center of (c) correspond to the bot- 
tom of the gray shading with an IPR value 7 = 0, whereas 
dark gray (red) patches inside the circular region in (c) and a 
dark gray dot well separated from bulk lying at extreme right 
correspond to the top of the shading scale with an IPR value 
/ — 800. The light gray vertical regions toward the left side 
in (a) and (b) is not as delocalized (with an IPR value close 
to 500, as in the scale) as the light gray vertical region toward 
the right side of the same plot (with an IPR value varying by 
about 300 in the scale). Similarly, light gray points at the 
right corner of (c) are more delocalized than the part of the 
vertical and horizontal gray area of similar shade crossing at 
center. Darker gray patches in (c) are not as delocalized as 
the patches with the same shade in (a) and (b). 

and both of them are delocalized. For r — 0.6, which 
means that nmxiber of directed connections is higher than 
T — 0.8, these isolated eigenvalues still exist on the real 
axis but come closer to the rest bulk part. Further 
increase in the number of directed connections brings 
isolated eigenvalues even closer, as well their localiza- 
tion properties get affected. For t = 0.4, these isolated 
eigenvalues are no more extended, and eigenstates hav- 
ing largest absolute values from group (A) and (B) are 
most localized. As r is decreased further, these isolated 
eigenvalues merges with the bulk part of group (A). 



VII. ZEBRAFISH METABOLIC NETWORK 

This section compares the behavior shown by the 
model networks discussed so far with a realworld net- 
work. We use network for the zebrafish metabolic system. 
The network is constructed by using the KEGG data base 
[3^ . Nodes in the network are reactions. If the product 
of one reaction (say A) is the substrate of another re- 
action (say B) then there is a directed link from A to 
B. There are also some directed links that correspond to 
the irreversible biochemical reactions. The largest con- 
nected component of the network has N — 1593 nodes 
and E = 9018 edges. The number of undirected edges 
is 6770 and the number of directed edges is 2248. The 



value of r can be calculated as r = 6670/9018 ~ 0.7 
The average degree of the network is E/N ^ 6. Fig- 
ure ^ plots the spectra of this network. In order to 
compare the results with a corresponding random net- 
work, we plot Figure ([9] c) spectra of a random networks 
ensemble with the same size and number of directed and 
undirected edges as in the zebrafish network. 

We notice that the spectrum of the zebrafish network 
has few drastic differences from those of model networks 
discussed in the preceding sections. The main observa- 
tion is the clustering of eigenvalues in an elliptical disk 
along the imaginary axis. Two extremal eigenvalues lie 
on the real axis, and they are far from the rest of the 
bulk eigenvalues. The corresponding random network 
also has the largest eigenvalue situated at the real axis 
with a value close to the average degree of the network. 
This eigenvalue is situated far from the rest of the bulk 
eigenvalues centered around zero. This is a typical spec- 
tral behavior of directed random networks, as discussed 
in the Sec. |IVl 

Furthermore, the localization properties of eigenstates 
of the real networks studied show very different phenom- 
ena than what observed for the corresponding model net- 
works. In an elliptical disk there is one layer of eigen- 
states that are very much localized, and the second layer 
in the disk shows a full range of IPR values between ex- 
treme localization and extreme delocalization. There are 
many eigenvalues that lie on the real axis, which can be 
explained by the relation between asymmetry in the ma- 
trix and the occurrence of complex eigenvalues discussed 
in Sec. VI (Fig. [8]). However the localization properties 
of these eigenstates corresponding to the real axis are dif- 
ferent than the corresponding model network and other 
earlier model networks we investigated in the preceding 
sections. The main difference is that all eigenstates lie 
toward the bottom on the localization scale. The isolated 
eigenvalues for model networks so far are observed with 
the most extended state, whereas for the zebrafish net- 
work, the isolated eigenvalue lies toward bottom of the 
localization value. 



VIII. CONCLUSION AND DISCUSSIONS 

We analyzed the spectral properties of various model 
directed networks and a realworld network. The network 
edges take value 1 and —1 depending upon whether it is 
starting from an excitatory node or from an inhibitory 
node. If all nodes are excitatory then the correspond- 
ing network is symmetric (r = 1). Directionality is in- 
troduced by making some nodes inhibitory, and conse- 
quently the corresponding edges take value —1. An equal 
expected value of inhibitory and excitatory nodes gives 
rise to a completely uncorrelated network {t — 0). 

Asymmetric networks (r = 0) with a community struc- 
ture modeled by a random sub-network show circular 
spectra similar to random uncorrelated matrices. Spec- 
tra for networks with entries 1,-1 follow Girko's circular 
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distribution with radius a/p]V(1~^~p). For a very well 
distinguished community structure Pr the whole 

spectrum is localized except for some of eigenstates at 
the boundary. As pr increases and the network deviates 
from a community structure there is a sudden change in 
the localization property for a very small value of defor- 
mation from the perfect community structure. The value 
as less as ~ 0.1, which corresponds to only 10% inter- 
community connections (Fig. [1]), gives rise to a delocal- 
ized spectrum. In order to understand this delocalization 
transition we track the spectra for very fine rewiring be- 
tween communities (Fig. [T|). For exactly one rewiring 
most of the eigenvalues lying on the real axis are delo- 
calized. This suggests that there are eigenfunctions con- 
fined within their respective random communities. With 
the increase in rewiring between communities some of 
the eigenstates lying at the boundary of the circular re- 
gion start delocalizing. This suggests that some of the 
eigenstates now span the whole network through rewired 
connections that connect two sub-networks. With the 
further deformation from the community structure, the 
next level of delocalization starts occurring through a 
circular ring around center, which gradually spans the 
whole spectrum with an increase in the rewiring. 

The spectra of random networks, where the probability 
for a node being inhibitory or excitatory is equal, show 
a circular distribution with a radius ^ypN{l — p). Based 
on the IPR values, the spectrum can be divided into two 
parts, (A) one part consisting of eigenstates at the real 
axis and at four corners with large absolute eigenvalues 
that arc localized, and (B) another part consisting of the 
bulk middle part of the spectrum, which is less localized. 
As the connection probability increases part (B) starts 
dominating the spectra, except for a few localized eigen- 
values, which remain very well separated from the bulk 
part even for very large connection probabilities. 

In order to understand the mechanism for localization 
we tracked the spectra as the network is rewired from a 
completely symmetric structure (r = 1) to a completely 
asymmetric one (r = 0). For symmetric networks spectra 
lie on the real axis with exactly one eigenvalue separated 
from the bulk. As connections are made directed by mak- 
ing some nodes inhibitory, some of the eigenstates start 
occurring in complex conjugate pairs. The eigenvalue dis- 
tribution along with the IPR value show a rich pattern. 
Overall, the spectra gradually become more delocalized 
as the number of directed connections increases. 

Furthermore we studied spectra of matrices with non- 
negative entries. This kind of matrices was previously 
investigated in details; we also showed that the spec- 
tral distribution follows Girko's generalized law. For 
T = it shows circular distribution, and as the corre- 
lation < AijAji > increases spectra become elliptical. 
The localization properties of whole spectra remain sim- 
ilar except for a sudden localization for very large values 
of r. 

Spectra of networks with a scalefree community struc- 
ture do not show a localization pattern that is much dif- 



ferent from spectra of networks with the random commu- 
nity structures. For two separated communities with one 
connecti on between them, the spectra are circular with 
radius \/pN{\ — p) and the boundary of the circular re- 
gion is scattered. As more connections are rewired lead- 
ing to stronger coupling between the communities, the 
spectral distribution remains same but the localization 
behavior of the eigenstates changes. The main difference 
from the spectra of networks in section Sec.|ll]is that they 
are less delocalized which is obvious because any network 
with a scalefree community structure always has more 
structure than the corresponding networks with a ran- 
dom community structure. The crucial point is that the 
spectral distribution is very robust with the changes in 
the structure of the network and depends strongly on the 
directionality measure t, whereas the localization prop- 
erty depends on both the structure of the networks and 
the directionality of connections, i.e., the value of r. 

We encounter interesting features while tracking the 
spectra as the network with a community structure goes 
from a directed to a symmetric one. For a random 
network with two random communities, the spectrum 
shows exactly two eigenvalues well separated from the 
bulk (Fig. [1]). These eigenstates correspond to the two 
largest absolute values and are delocalized. The largest 
one is completely delocalized as was the case for a random 
network without a community structure, but the second 
largest one that is well separated from the bulk provides 
a distinguished signature of two communities in the net- 
work. As the number of directed connections increases 
the spectrum can again be divided into two parts, (A) 
the part of the spectrum with real eigenvalues and (B) 
the part of the spectrum with complex conjugate pairs. 

To the end we analyze a realworld spectrum. Spec- 
tral properties of the metabolic network of zebrafish show 
much richer features than the model networks considered 
in the present paper. One very striking difference is the 
localization property of the eigenstate with the largest 
absolute value. This eigenstate is localized in contrast to 
those of the model networks. Another observation is the 
clustering of eigenvalues in an elliptical disk along the 
imaginary axis. 

In the paper we have investigated spectra and localiza- 
tion properties of directed networks with binary entries. 
The networks with inhibitory and excitatory nodes have 
much richer spectra than the networks with only excita- 
tory nodes. The bulk of the spectra for completely asym- 
metric networks follow Girko's law, but as probability of 
directed connections is reduced the spectra show very 
different patterns depending upon the network structure 
and the ratio of inhibitory and excitatory nodes. The 
realworld network that we studied here shows a very dif- 
ferent spectral pattern from any of the model networks we 
have studied. Though directed networks span a variety of 
complex systems, the research for directed networks lead- 
ing to complex eigenvalues is limited. The results pre- 
sented in the present paper provide a useful platform to 
understand the structural pattern in directed networks. 
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and can be used to investigate further the dynamical be- for providing zebrafish metabohc data from KEGG [3! 
havior of nodes relevant to a variety of problems ranging 
from physics to sociology. 



IX. ACKNOWLEDGMENT 



We acknowledge Dr. Changsong Zhou for useful dis- 
cussions at the initial phase of the work, and Zhang Xun 



[1] M. L. Mehta, Random Matrices, 3rd ed. (Elsevier Aca- 
demic Press, Amsterdam, 2004). 

[2] T. Guhr, A. Muller-Groeling and H. A. Weidenmuller, 
Phys. Rep. 299, 189 (1998). 

[3] R. Grobe, F. Hakke and H.-J. Sommers, Phys. Rev. Lett. 
61 1899 (1988). 

[4] G. Palla and G. Vattay, New J. Phys. 8 307 (2006). 

[5] J. N. Bandyopadhyay and S. Jalan, Phys. Rev. E 76, 
026109 (2007). 

[6] S. Jalan and J. N. Bandyopadhyay, Phys. Rev. E 76, 
046107 (2007). 

[7] M. A. M. de Aguiar and Y. Bar- Yam, Phys. Rev. E 71, 
016106 (2005); A. N. Samukhin, S. N. Dorogovtsev, and 
J. F. F. Mendes ihid 77, 036115 (2008). 

[8] K. -I. Goh, B. Kahng, and D. Kim, Phys. Rev. E 64, 
051903 (2001); S. N. Dorogovtsev, A. V. Goltsev, J. F. 

F. Mendes and A. N. Samukhin, ihid 68, 046109 (2003); 
E. Estrada, Europhys. Lett. 73 (4) 649 (2006); D. Kim 
and B. Kahng, Chaos 17 026115 (2007); H. Yang, C. Yin, 

G. Zhu and B. Li, Phys. Rev. E 77 045101(R) (2008); G. 
Zhu, H. Yang, C. Yin and B. Li, Phys. Rev. E 77, 066113 
(2008); G. Bianconi, arXiv:0804.1744. 

[9] S. Jalan, N. Solymosi, G. Vattay, B. Li, Phys. Rev. E 81, 
046118 (2010). 

[10] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[11] D. J. Watts and S. H. Strogatz, Nature 440, 393 (1998). 

[12] S. Jalan, Phys. Rev. E 80, 046101 (2009). 

[13] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002) and references therein; S. Boccaletti, V. Latora, 
Y. Moreno, M. Chavez, D.-U. Hwang, Phys. Rep. 424, 
175 (2006). 

[14] M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. 

USA 99, 7821 (2002); M. E. J. Newman, Social Networks 

27, 39 (2005); M. E. J. Newman, Proc. Natl. Acad. Sci. 

USA 103, 8577 (2006); M. J. Krawczyk, Phys. Rev. E 

77 065701 (R) (2008). 
[15] E. Ravsaz et al, Science 297, 1551 (2002); R. Guimera 

and L. A. N. Amaral, Nature 433, 895 (2005). 



[16] G. Palla, I. Derenyi, I. Farkas and T. Vicsek, Nature 435, 

814 (2005); M. E. J. Newman, Phys. Rev. E 70, 056131 

(2004); A. Arenas, A. Fernandez and S. Gomez, New J. 

Phys. 10 053039 (2008). 
[17] V. L. Girko, Theory Probab. Appl. 29 694 (1984). 
[18] H. J. Sommers, A. Crisanti, H. Sompolinsky, and Y. 

Stein, Phys. Rev. Lett. 60, 1895 (1988). 
[19] O. Giraud, B. Georgeot and D. L. Shepelyansky, Phys. 

Rev. E 80, 026107 (2009) 
[20] B. Georgeot, O. Giraud, and D. L. Shepelyansky Phys. 

Rev. E 81, 056109 (2010). 
[21] T. Rogers and I. P. Castillo, Phys. Rev. E 79 012101 

(2009). 

[22] J. J. Seidel, Lin. Alg. Appl. 1 281 (1968). 

[23] J. J. Seidel, A survey of two-graphs, 481 Proc. Internat. 
Coll. Teorie Combinatone (Roma 1973), Accad. Naz. 
Lincei, Roma, (1977); J. J. Seidel and D. E. Taylor, Two- 
graphs: a second survey, 689 Algebraic methods m graph 
theory, Vol. II, Conf. Szeged 1978, CoUoq. Math. Janos 
Bolyai 25 (1981). 

[24] Course notes by A. E. Brouwer and W. L. Haemers, avail- 
able at http://www.win.tue.nl/ aeb/math/ipm.pdf 

[25] J. G. Martin and E. R. Canfield, SIAM J. Sci. Comput. 
31 5, 3529 (2009). 

[26] K. Rajan and L. F. Abbott , Phys. Rev. Lett. 97, 188104 
(2006). 

[27] P. Erdos and A. Renyi, Publ. Math. Inst. Hungar. Acad. 

Sci. 5, 17 (1960). 
[28] O. Hul, P. Seba and L Sirko, Phys. Rev. E 79 066204 

(2009) 

[29] B. BoUobas, Random Graphs, (Academic Press, London 
[30] F. Juhasz, Discrete Math. 41 161 (1982). 
[31] D. W. Lang, Physical Review 135 Number 4B, B1082 
(1964). 

[32] Kyoto Encyclopedia of Genes and Genomes (KEGG); 
http://www.genome.jp/kegg/ 



